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On a Certain Completely Integrable System of Linear 
Partial Differential Equations. 

By E. J. WlLCZYNSKI. 



Introduction. 



From the point of view of projective differential geometry, developables 
and ruled surfaces require a special discussion, on account of their many 
exceptional properties. These exceptional properties manifest themselves 
analytically by the fact that, for these surfaces, the differential invariants of 
lowest order are equal to zero. If a surface is neither a developable nor a 
ruled surface, its absolute invariants of lowest order are those which have 
been denoted by the author by I and J. In a certain sense, then, those surfaces 
for which I and J are identically equal to zero, are to be regarded as consti- 
tuting an especially simple and fundamental class. 

I have shown recently,* that these surfaces are the integral surfaces of 
the completely integrable system of partial differential equations 



(8) 






where c , c x and c 2 are constants, and I have studied in detail the case c = 0. 
The system (8) is then capable of integration by means of elementary functions, 
and the corresponding surfaces, for which, of course, I=J=0, have the further 
property of being invariant under a continuous group of projective transfor- 
mations. If c ^=0, no such group exists, which leaves the surface invariant, 
and the problem of integrating (8) can not be solved in such a simple manner. 
The directrix curvesf of any integral surface of (S) form a conjugate net, 
whether c vanishes or not. Upon this remark may be based a transformation 
theory of system (8) which we hope to discuss on some future occasion. The 

* "On a Certain Class of Self -projective Surfaces," Transactions of the American Mathematical 
Society, October, 1913. 

f For the notion of directrix curves, cf. a paper of mine in the Transactions of the American 
Mathematical Society, Vol. IX (1008), pp. 114-120. 
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present paper, however, is merely devoted to the integration of system (8) in 
the case c zf= 0. 

In order to accomplish this integration, we find it necessary to solve the 
following problem: To integrate a certain completely integrable system of 
non-homogeneous linear partial differential equations, when the solutions of 
the corresponding homogeneous system are known. Since I have been unable 
to find any discussion of this question in the literature, the first two sections 
of this paper are devoted to a solution of this problem in so far as it is needed 
for what follows. It then becomes possible to prove the existence of various 
kinds of solutions of system (8). The most interesting of these can be written 
as an exponential function, multiplied by a series of positive integral powers 
of c , the coefficients of this power- series being polynomials in u and v. 
A number of methods will be given for calculating these polynomials and for 
studying spme of their properties, but I have been able to find a perfectly 
explicit formula for only some of them. 

The Solutions of (8) are found to satisfy integral equations of a peculiar 
kind, involving integrals of exact differentials. A more detailed discussion of 
such integral equations will be given in a separate paper. 

§ 1. Integrability Conditions for a Non-Homogeneous System of Linear 
Partial Differential Equations of the Second Order. 

Let us consider the non-homogeneous system of partial differential equations 



where 



(L) { 



*/«« + 2ay u + 2by v + cy + d = 0, 
y vv + 2 a'y u + 2 b'y v + c'y + d'=0, 



_dy dy d 2 y d z y 

Vu ~du' Pv ~dv' VvM ~du^ Vvv ~W 



and where a,b, . . . . , d' are analytic functions of u and v. By differentiating 
these equations, we obtain, for the third-order derivatives of y, the following 
unique expressions: 

Vuuu — PlVm + PiVu + VsVv + PiV + Pn , 

Vuuv — Q-xVuv + &y« + & Vv + liy + & t 

Vuw = r t y uv -+ r 2 y n + r % y v + riy + r h , 

Vvw = s i Vuv + s 2 y u + s 3 y v + s t y + s 5 , 
where 

g x = — 2a, q 2 = 4:a'b — 2a v , q 3 = 4=bb' — 2b v — c, q i = 2bc' — c v , q 5 = 2bd' — d v , 
r x = — 2b', r 2 = 4caa' — 2a' u — c', r 3 = 4:a'b — 2b' u , r i = 2a'c — c' u , r 5 = 2a'd — d'„, 

and where the values of p t and s t are immaterial for our present purpose. 
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If we form the fourth-order derivatives of y, the expressions for y muu and y, 
are unique, and the two expressions found for y unm , viz., —Mum and x uuv , 



u««v, *—, dv — du 



will be identical without imposing any further condition, since each of them 
is the result of differentiating the first equation of (L) once with respect to u 
and once with respect to v. Similarly, the equation 

® y«w & y vw 



dv du 

will be satisfied without imposing any conditions on the coefficients of (L). 

d '14 O 14 

The equality of the two expressions J? uuv and ;p OT , however, gives rise to 
such conditions. In fact we find 

^§T = (If + ^ + & ) y - + (If + «i r «- 2a '2«) y* 

+ (If +2i r «- 26 '««+«0 y.+ (|| +3i»-4-C&) 2/+ |f +air 5 -d / 2„ 

+ (l¥ + riS3_2&r2 )^ + (IS + ri ^ _cr2 ) 2/ + IS + rie5_ir2 

In order that system (L) may have four linearly independent solutions, 
the corresponding coefficients of y uv , y u , y v , y and y° in these two expressions 
must be equal.* We obtain in this way the following five integrability conditions : 



(I) 



a v — K = 0, 

aL + c' u — 2 a' a u — 2 a a' u — (a vv + 2 b' a v — 2 6 a' v — 4 a'. & J = 0, 

&:„ + 2afc: — 2a' &„ — 4 &«,; — (&„ + c, — 2 6 &; — 2 &'&,) = 0, 

c^ — 4cai — 2a'c, + 2a< — (c„ — 4c'6, — 2 6c; + 26'c,) = 0, 
J' um + 2adl—2a'd i , — 4 ! a' u d — c'd — (d„ + 2b'd v — 2bd , v — 4 : b v d'--cd') = 0. 



Thus, the integrability conditions for the non-homogeneous system (L) 
differ from those for the corresponding homogeneous system only by the 
presence of the last of these five conditions, t 

* Otherwise y would also satisfy an equation of the form 

yuv + 2 a" yu + 2 o" y v + c" y + d" = 0. 
The system obtained by adding this equation to (L) has at most three linearly independent solutions. 

f E. J. Wilczynski, " Projective Differential Geometry of Curved Surfaces " (First Memoir) , 
Transactions of the American Mathematical Society, Vol. VIII (1907), p. 246. 
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§ 2. Integration by Quadratures of the N on-Homogeneous System in Terms 
of the Solutions of the Corresponding Homogeneous System. 

If the conditions (I) are satisfied, (L) is a completely integrable system, 
and the corresponding homogeneous system 

Vuu + 2 a Vu + 2 b y v + c y =■ 0, 
I y vv + 2a'y n + 2 6'y, + c'y = 

will also be completely integrable. Therefore, (H) will possess precisely four 
linearly independent solutions. Let y (1) , y (V >, y (S) , -y (i) be four such linearly 
independent solutions of (H) , so that 

c x yV + c 2 y<v + c s y™ + c t y& 

will be its general solution if c x , . . . . , c 4 are arbitrary constants. Let y be 
the general, and Y any particular solution of (L) Then y — Y will be a solu- 
tion of (-BT)j so that 

y = Y + c l2/ w + c 2 ?/< 2 > + c 3 */< 3 > + c 4 y<*> 

Consequently, if we assume the solutions of (H) as known, it suffices to find 
any particular solution of (L) in order to integrate (L) completely. 

We use the method of variation of constants to find such a particular 
solution of (L). Let us write Y in the form 

Y = Xl (k) y (k \ (1) 

where Z (1) , . . . ., l (i) are functions, yet to be determined, of u and v. We shall 
have 

if we impose upon Z< 1} , . . . . , Z (4) the conditions 

2Z<*V* ) =2^*V* ) =0, 
or 

2y<*>d*<*> = 0. (2) 

We then find 

Y uu = 2 (/<*> yS> + Z<*> y<*>) , I,, = 2 (Z<*> y<*> + «*> 3/? >) , 

so that Y will be a solution of (L), if, and only if, 

2 *<» *,<*> = - d, 2 J<*> #<*> = - d'. (3 ) 

We may write 

Y = 2 l <k) u (k) 

if we impose the further conditions 

2*«3tf> = 0, 2l*>yW = 0. (4) 



(6) 
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But equations (3) and (4) are equivalent to the following: 

2y™dl™ = —d-du, 2y™dl™ = —d'-dv, (5) 

so that we have in (2) and (5) a system of three equations for the four 
differentials dV k) . We add to these equations a fourth, viz., 

2y&> d*<*> = — {Adu + Bdv) 
where we shall determine the (as yet unknown) functions A and B of u and v 
in such a way that the values of dl^ obtained from these four equations 
shall be exact differentials. 

Let us write these equations as follows : 

2y£>dZ<*> = — (Adu + Bdv), 
Xy™dl™ = —d-du, 
2yi k *>diw = —d'-dv, 
Xy m dl<^ = 0, 

and denote the determinant of the left members by A, so that 

A= | y„, y u , y v , y j 
Then A is not identically equal to zero. For, if it were, (H) would have at 
most three linearly independent solutions, and y( l) , . . . . , y w could not be 
linearly independent. 
"We may write 

A = 2 y™ *,<*> = 2yi k) p<*> = 2 yi h) v™ = 2 */<*> p<*\ 
where ^ (ft) , ^ (ifc) , etc., are the third-order minors of A, and may be written 

as follows: 

^d) = + \y u , y v , y\ f ^d) = — \y„ v ,y v , y\, 

v w = + \y uv , y u , y\, p (1) = — \y uv , y u , y v \, 
etc., the notation adopted being self-explanatory. 
The solution of (6) for dl m gives 

Adl^ = — X<*> (Adu + Bdv) — p<*> d du — v^d' dv 
But we find 

W = \y*u, Vv, y\ + \y u , y™, y\ + \y u , Vv, Vu\ 

= \ — '*ay u — 2by v — cy,y v ,y\-+ \y u , y m ,y\ =-2a^)-r« 
and similarly, 

a4*>.-= — 2aa<*> — *<*>, a,« = — 2 6'^w— ,u<*>, (7) 

whence 

p<*> = — a<*> — 2 b' K*\ v w = -^>-2a *,<*> 

Consequently we may write 

Adl™=L™d« + M™dv, (A; = 1,2,3,4), (8) 
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::} 



(9) 



if we put 

L<»> = d •*,<*> + (2 6' d — A)M*>> 
M w = A , . X (*) + ( 2 a # _ B) a/*>. 

In order that the expressions for dZ w may be exact differentials, we 
must have 

aA A / 3«A A / 
If we make use of the equations (7) and the equations 



3A 



= — 4 a A, 



3A 



= —46' A, 



3 m ■*■<*—» g ^ 

we find that these integrability conditions give rise to a system of four equa+ions 
of the form 

(B — d' u + 2 a' d) *,« — (4 — d„ + 2 & d') Jl<*> 

+ [— A v + B u + 2b'd v — 2ad u — d(2a' u + C — 8 a a') 

+ d'(26„ + c — 8 6 6') — 4 6'-4+4aB]X<*>=0, (10) 

(A; = 1,2, 3,4). 

If we interpret 2/ (1) , . . . . , ?/ (4) as the homogeneous coordinates of a point P & 
in space, the point P y describes a surface as u and v pass through the values 
of their respective ranges. Moreover, this surface can not be a developable, 
for a system of equations of the form (H ) . f But the quantities yl ( *> are the 
homogeneous coordinates of the plane, tangent to this surface at P y , and 
therefore can not satisfy a system of equations of the form (10) with non- 
vanishing coefficients. For, if they did, their ratios would be functions of a 
single variable; i e., the surface would be a developable. Consequently, the 
coefficients of (10) must be equal to zero. If we equate to zero the first two, 

we find 

A = d v — 2 6 d', B = d' u — 2a'd, 

and these values of A and B also cause the remaining coefficient to vanish, as 
a consequence of the last of the integrability conditions (I). 

We have obtained the following theorem : 

Let y a) , . . . . , 2/ (4) be four linearly independent solutions of the homo- 
geneous system (H) , and let 



(11) 



A = 


U W u (2) 
Vu y% 

tf> yT 

y m y m 


i)wo 

yi 3) 
y (3) 


i/uv 

yi iy 
yV 




* hoc. oit., 


p. 258. 






t Ibid., § 2. 
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so that /l (fc) is the co-factor of y$ in this determinant, which, moreover, is not 
identically equal to zero. If we put 

L« = d X<*> + (2 6 d' + 2 V d — d v ) V*\ 
M <*> = d'K k) + (2 a d' + 2 a' d — d' u ) *<*>, 
the quantities 



] (ft = 1,2, 3, 4), (12) 



dl^= 1 (LMdu + M^dv) (13) 



wiW fee e#ac£ differentials, and the function 

Y = Xl w y w =iy (k) I - (L^du + M^dv), (14) 

k=l Jt=l / &■ 

where the lower limits of the integrals are arbitrary constants, will fee a solution 
of the non-homogeneous system (L). 

If we choose the four lower limits equal, so that 

a k = a, b k = 6, (lc = 1, 2, 3, 4) , 

we obtain that particular solution of (L) which corresponds to the initial 

conditions 

Y = Y U = Y V = Y, w = f or u = a, v = 6. 

The quantities X (fc) satisfy a system of partial differential equations of the 
same type as (JET). If we put 

r w =^, (4 = 1,2,3,4), 
J (1) , , T (4) are solutions of the adjoined system of (jff), viz.: 

r r tm + 2ar u -2&r„+(c + 2&„-4&&')r = o, 

*" ' I r w -2aT.+ 26T.+ (c'+2a:-4aa')r = 0.' 
In terms of these quantities we find 

L m _ va [d r<*> + (2 fed' — dj I«] , 

M(fc ) = va [dT« + (2 a'd — dl) 7< fc >] , 
and these formulae are sometimes more convenient than (12). In all appli- 
cations of these formulae it is good to remember that, since 

du ov 



(12') 



* Loo. cit., p. 259. It should be remarked that system (/?') is not in general equivalent to the 
system obtained from (ff) by taking the Riemannian adjoint of each of its equations. There are two 
distinct methods of generalizing the Lagrange adjoint of an ordinary linear differential equation. The 
above method, based on the principle of duality, may be called the geometric method, as distinguished 
from the analytic method of Riemann. 
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we shall have 

A = C , e~ 4 *', C = const., (15) 

where p is determined by the conditions 

^£ - a ?P. - b' 

which are consistent, according to the first of the integrability conditions (I) . 
The constant C may, of course, be reduced to unity, by a proper choice of the 
fundamental system y^\ . . . . , t/ (4) 

§ 3. Integration of System (S) by Series Proceeding according to Positive 

Integral Powers of c 

If the coefficients of an ordinary linear homogeneous differential equation 

d n y , t> d n ~ x y . . D n 

are integral rational functions of a parameter (i, and if, for an initial value of 
the independent variable x = x independent of fi, initial values be prescribed 
for y, y', . . . ., t/(" _1) which are either independent of (i or integral (rational 
or transcendental) functions of (i, the corresponding solution of the differential 
equation will itself be an integral (rational or transcendental) function of p. 
This theorem, due to Poincare, Gunther and Horn, may be applied to our 
system, in so far as we may by simple processes find from (S) an ordinary 
differential equation of the fourth order for y as a function of u, the coefficients 
of this equation being rational integral functions of c and v, as well as of u 
A similar equation of course exists for y as function of v. The actual deter- 
mination of the solution might be accomplished on this basis, but manifestly 
the lack of symmetry in the treatment accorded to the two independent variables 
would be a disadvantage. We therefore prefer to start anew, without making 
any use of the above theorem except in so far as it assures us a priori of the 
existence of solutions of the form 

y r = y? ) + y?'c« + y? ) cl+ ...., (16) 

and even this is unnecessary, as we shall establish the existence of such solutions 
directly a posteriori. 

Let us then assume for y in system (S) a development of form (16) We 
find at once the following conditions : 

* Poincare, Acta Mathematica, Vol. TV (1884), p. 212. Gunther, Crelle's Journal, Vol. CVII (1889) 
p. 312. Horn, Mathematisohe Annalen, Vol. LII (1898), p. 343. 
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a^> V <o)_ 



(17a) 



and 



3 2 « (4) 3 « (ifc > 

° ff f 4- 2 ^ r 4- c «<*> 4- v i/* -0 — 
dv 2 du ~ 



(fc = l,2,3,....), (17b) 



which, together with the proper convergence conditions, are necessary and 
sufficient for the existence of a solution of (8) of the form (16). 

According to (17a) we may take for yf^> the exponential e arW+PrV , provided 
the constants a r and /? r satisfy the equations 

a» + 20 + Cl = O, "I 
£« + 2a + c, = 0. J 
In general, there are four linearly independent functions of this kind, 

j,fo) ==e . i „ + /».. > (4 = 1,2,3,4), (19) 

and we shall, for the present, confine our discussion to this case. As we shall 
see at the end of this paragraph, this limitation does not essentially affect the 
generality of our argument and, moreover, may be easily removed afterward 
The solution yi 0) = e arW+p,v of (17a) is characterized by the initial conditions 

r* (o) Q (0) CJ2 (0) 

■w (0) = g«r«+/3r6 °Vr __ a e a r a + p r v g Vr _ O e a r a + r b ° Vr _ a O e a r a + p r b 

' ' du r ' dv Hr ' 3w3v rHr 

for w = a, v = 6. We shall determine that solution of (S) which corresponds 
to these same initial conditions. If this solution is expressible by a series of 
form (16), we must satisfy the initial conditions 

3« (fc) dii w 3 2 w (fc > 
du dv duov — 

for u = a, v = b, while 

We proceed to calculate z/< ft) by applying the method of § 2 to (17b) In 
the notation of § 2 we have 

a = 0, 6 = 1, c = c l , d = uy$ k ~ 1) , 
a' — l, 6' = 0, c' = c t , d' = vy^- 1) 
Of the five integrability conditions (J), the first four are obviously satisfied. 
The left member of the last reduces to 

* A complete discussion of the various cases is to be found in my paper " On a Certain Class of 
Self-Projective Surfaces," to which I hare already referred. 
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11 ( ° y r A- 9 y r 4- r «<*-!) 1 u( ^ r -X- 2 V r I- r 7/ (fc_1) 1 

and this is equal to zero on account of (17b), if k>2, and on account of (17a), 
if k = 1. 

The homogeneous system, which corresponds to (17b), is (17a), and its 
solutions are determined by (18) and (19) Therefore, a x , . . . ., a 4 are the 
four roots of the equation 



a* + 2 c x a 2 + 8 a + c\ + 4 c 2 



0, 



(21) 



so that 

ol x + <x 2 + .<x 3 + a. = & + & + ft + & = 

Consequently we find 

yf stf» yf> 2/i 0) - 3 
and 

«lA a 2^2 « 3 Aj «4& 



A A A A 

1111 



a? 


„ 3 

a 2 


™ 3 


«J 


«i 


a 2 


a 3 


a 4 


a? 


«S 


a| 


«! 


1 


1 


1 


1 



(22) 



or 



4 A = — (a 2 — a 2 ) (a x — a 3 ) (a x — a 4 ) (a 2 — a 8 ) (a 2 — a 4 ) (a 3 — a 4 ) 

We then find, for the co-factor of ~ &X in A, the value 

dudv 



W-*> = 



2yr* 



(* = 1,2,3,4), 



(23) 



where 

h x — (a 2 — a 3 ) (a 2 — a 4 ) (a 3 — a 4 ) , & 3 = (a 4 — aj (a 4 — a 2 ) (a x — a 2 ) , 
& 2 = — (03— a 4 ) (a 3 — aj (a 4 — <Xj) , & 4 = — (04— a,,) (a,— a 3 ) (o^— a 3 ) 

If we substitute these values in (12), we find 

AfJfJ = ** e — «— A" f"(2 u—a t v) yf'^ — v dy ^ ~\ , 

where we have introduced a slight change of notation as compared with (12), 
the upper index k of (12) having been replaced by the lower index i, while the 
significance of the two new indices r and k, which appear in the last two 
equations, is obvious. 

From the above values of A and k x we obtain the value of their quotient, 

4A 



h 



— — (a,— a 2 ) (aj— a 3 ) (04— a 4 ) 
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But, on account of (21), we have the relations 

2(^ = 0, Xa i a f =2c l , Xa i a i a k = — 8, a 2 a 2 a 3 a 4 = c\ -f- 4c 2 , 



so that 

4A 



^ = -4(a? + c 1 a 1 + 2) = 8(a x & — 1), 

proving the first of the four equations 

| = 2(a,fr-l), (4 = 1,2,3,4) (24) 

Cleaiiy, a^ — 1 will be different from zero, since a s is a simple root of (21) 
From (14) we find, finally, 

y m = 2e a "' +ft »/ , e ~ a ^~^!, I I ( — ft w + 2 v ) y? ' 5 > - u --Mf : - 7 ~- \ d u 

a, b 



4 / p-aiu-p+v re , 3«(* _1 >1 

*«~ + *7 < M r i ) [{'-ft' +i "^-irj 

*^ a, b 

+ U2u — * i v)y?-» — v dy £~" '\dv\ (25) 

The lower limits of the four integrals which appear in this formula have been 
chosen equal to (a, b) so as to satisfy the conditions (20) That the expression 
(25) for yf^ actually satisfies these conditions may be seen as follows. We 
may write (25) in the form 

y™ = h<Ry?>, (25a) 



where the quantities V*\ have been determined, in accordance with the method 
of § 2, in such a way that 

a Vr y 7(fc) ° Vi a Vr y 7(fc) Oft ° ffr y 7(fc) ° Vi 

du ~ i t 1 r ' i du ' dv ~ i = 1 r ' i dv ' a^a^ - ^ r > r dudv 

Consequently, all of these expressions will actually reduce to zero for u = a, 
v = b 

In order to establish the convergence of the expansion (16) for y r , it will 
be convenient to write 



K i (u, v; p, q) = 



8(a,&-l) ' 



(» = 1,2,3,4), (26) 
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Then we shall have 



Ml 

du 



/u, v 

4 r „ ^ a w »-» a,,(»-i) 



a_J/L fc) _ 

dv 



.2 ft *, [*<*-» <Z £, - £§^- d(p>)- °%— d ( fl «) j 

«^ a, b 



(27) 



3m 



A simple induction proof, involving the use of equation (25), shows that, for 
every value of Jc, yi k) and all of its partial derivatives are continuous (and even 
analytic) functions of u and v, for all finite . values of these variables. The 
path of integration from (a, b) to (u, v) may be chosen arbitrarily in the finite 
portion of the plane, since all of the differentials involved in (25) are exact. 
We may, moreover, without any essential restriction of generality, assume 
a = h = 0. 

Let N be a positive number not less than the largest value assumed by 
| e a rP + ^ r «| on fljg path of integration. We shall then have 



M<N, 



dp 



£\*r\#, 



-3g 



£\Pr\N- 



For any assigned finite values of u and v, there will exist a positive con- 
stant M 4 , such that, for all values of p and q along the path of integration, 

\K t (<u,v; p,q)\ <M t . 
Let 

M= 2M f , 



i=l 



and denote by a and /? two positive numbers, such that 

|a,|<a, |ft|<0, (i = 1,2, 3,4), 

and consequently, 



2 K t (u, v;p,q) 



<M, 



2a^ 



<al, 






<0M. 



Introduce three positive functions by putting 

P= C P ' 9 \d(p*)\, Q= C q \d{q>)\, 8= C 9 d^\p\^4,\p\\q\+o.\q^).. 

•SO, " 0, '-'O, 
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and we shall have 

\dS t \<dS, (t = l,2,3,4) 

If now we put h = 1 in (27), we find 

J. 11, V i , 

2, Mt\N dS + (3N dP + a.N dQ\ <N M T, 
0, i=l 

where 

and similarly, 
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O) 



du 



<aNMT, 



(i) 



dv 



<(3NMT, 



(i) 



dudv 



<a(3NMT. 



If we make use of these inequalities after putting & = 2 in (27), we find 



,T,f2 7^2 

|^)|<^^_L, etc., 



and, by induction, in genera], 



\y^\<N 



Mill 

Jcl 



dy ( r 



(fc) 



3w 



<aN 



2 
A!" ' 



(*o 



a.v< 

3 t> 






These inequalities assure us of the absolute and uniform convergence 
of our series for y r , and of its term-by-term differentiability with respect to 
u and v taken either separately or together. It is easy to show, further, that 



may also be obtained by repeated term-by-term differentiations 



Pi and Pi 

or ov z 

from y r as absolutely and uniformly convergent series. Consequently, we have 

obtained the following theorem : 

The series 

GO 

y r = 2 yj*>e$, 

/or which yf^ = e arU+erV , and whose general term is determined by the recursion 
formula (25), is absolutely and uniformly convergent for all finite values of 
u, v and c , and represents that solution of the completely integrable system 
(S) , whose initial values for u = a, v = b coincide with those of the exponential 
function e arU+firV 

We have computed the solutions y r under the assumption that the four 
roots a 1} . . . . , <x 4 of, (21) are distinct. But this assumption, as we have already 
stated, does not really restrict the generality of our method. For, by means 
of a linear transformation of the independent variables 

u — u-\-h, v = v-\-Jc, (28a) 

system (S) may be transformed into a new system of the same kind, whose 
constants have the values 
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^o^Co, d 1 — c 1 + hc Q , c 2 = c 2 + k c . (28b) 

If c =/bO, we may therefore transform (8) into another system of the same 
form whose constants c x and c 2 have any desired values, for instance c x — c 2 = 0, 
in which case the four roots of (21) will actually be distinct. On the other 
hand, it would not be difficult to develop formulae to take the place of (25) in 
the case of two or more equal roots. 

The linear transformation (28a) enables us to materially simplify the 
appearance of the differential equations (8) , but we have made no use of this 
transformation in developing our theory up to this point, for two reasons. In 
the first place, there would have resulted no essential simplification in any of 
our formulae sufficient to make up for the corresponding loss of symmetry. 
But, in the second place, the general formulae have a real significance which 
might otherwise have escaped our notice, and which gives us a real insight 
into the nature of the integral surface. 

In fact, the linear transformation (28a) may be interpreted as a change 
of position on the integral surface, instead of as a change of the parameters of 
reference. Now, if c is small, our formulae show that the surface whose 
parametric equations are 

y r ^e^ n+ ^ v , (r = l, 2, 3,4), (29) 

may be regarded as protectively equivalent to a first approximation of our 
integral surface in the neighborhood,, of the point for which u = a, v = b. 
If a r and ~R r are the values of a r and /? r which correspond to the constants 
Cj and c 2 , the surface 

y r — e 5 r « + A-», (r = 1, 2, 3,4), 

will, in the same way, be protectively equivalent to a first approximation of 

the integral surface in the vicinity of the point u = a -{- h, v = b + h. In other 

words, for small values of e , an integral surface 2 of (8) has associated with 

each of its points P, a self -projective surface of the form (29), whose exponents 

a r , @ r change with P, and which may be regarded as a first approximation of 2 

in the vicinity of P. 

One further remark is of importance for our later developments. From 

the nature of our convergence proof by dominant functions, it is clear that the 

series 

yP + ypc + ypcl+ .... 

will still be a uniformly and absolutely convergent series for all finite values 
of u, v and c , if the recursion formula (25) be altered by writing in place of 

* Of. foot-note on page 239. 
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{a, b), the common lower limit of all of the integrals, a different lower limit 
{a\ h) , bf^) for each of the integrals, provided that all of these limits are finite. 
Moreover, the function of u and v thus defined will still be a solution of (S), 
corresponding, however, to altered initial conditions. 

§ 4. Properties of the Coefficients of the Series Obtained in the Preceding 

Section. 

Let us make the generalization referred to at the end of the preceding 
section, and on that assumption calculate yp. If we make use of the abbre- 
viation involved in writing (25) in the form (25a), we find 

4(a r /5 r — 1) L J (a,o>, &,<•>) 

9T7TR — TT7« 7TTTR Tpr S («r~ a,)w + (^,— A)v — IN ,i=t=r. 

Then 

,/d _ ■,«,«+/!,«[ - /g y w 2 + 2«i;— oc r ?; 2 * (a,-— optt +(&—&) ^—t "I 
^ " L 4(a r/ 3 r -l) ' ^^2(^^-1) (a ,-«,) (fr-ft) J 

4 <,(a,-a,)ai< 1 >GS r -/SOW , > 

- 2 ' e °'" +/M o/„/g iw„ „w/i /^ S («r-«i)«. (,) + (/» -A)6? } -1 

i=i 2,(a i p i -~±){a. T —a i )(p r —p i ) 

W^l) ' (30) 

where the stroke on the summation symbol indicates that the value * = r is to 
be omitted. It will obviously be possible to find finite numbers a$ l \ 6j 1) and 
a?\ b^ so as to cause the terms e aiW+PiV (i=f=r) to disappear, and to reduce to 
zero the constant term of the expression which is multiplied by e arU+l3rV If 
this be done, we have 

,/D _ p a T «+f) r v r —{3 r u 2 + 2uv—a r v* * (ar—aJit+iPr—PJv ~| 
Vt ~ e L 4(a r fr-l) ^ ZiZioLiP-1) {*-*>) (Pr-Pi)\ 

= 6 M+fc, P{"(« ) »), (30 a) 

so that yl X) reduces to an exponential factor multiplied by a polynomial Pp (u, v) 
of the second degree in u and v, whose constant term is equal to zero. 

We may generalize this result. For every value of k there exists a solution 
yf^ of (17b), expressible in the form 

2 k 

where .4$ is a homogeneous polynomial of degree j in w and v. In other words, 
yj k) is equal to e°r«+A-i> mu itiplied by a polynomial of degree 2k in w and v, 
without a constant term. 
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The theorem is true for k = l. We proceed to prove the general theorem 
by induction. Let us suppose, therefore, that 



and consequently 



du 



2fc 2 



^ = e a r u + Pr v X (<**?>*- + p r A?T*) . 



We shall find 



. («, v) 



a = f e (-r--o.+w,-ft)« [ |p<*-i) + p«-y +....+ pft-ft-i | d M 

+ 1 G&P +■■••• + Qfr&-i S ^ v] , (31) 

where 

P ^ = 4(«J,-1) [l-(/gr + ft^ + 2^1 AfrJl-u^fip], 



(fc-i) 



f)(fc-r 



4( 



«*ft-i)L {: 



3v 

i T>1 



]■ 



(32) 



l2u-{* r + a. i )v\A?jll-v'^ u 

(j = l,2>,3,....,2Jc-l), 
in the application of which formulae we must remember to put 

A^ = Afc$ = 0. 

Clearly, Pfj^ and Ql k ,i, l j are homogeneous polynomials in u and v, of degree j. 
From our general theory we know that the expression under the integral sign 
in (31) is an exact differential. We must, therefore, have 

(ft, - &) Pftft-i = (a r - a,) GftiU , 

\Pr Pi) r r,i,2k-2 ~\ p~T. — \ a r a i) Vr,i,2fc-2 ~i o— 



-1) 
,2fc-l 



dv 



<u 



■> p(fc-i) a/i(*-i) 

(&-&) Pft? + ^T 1 = (a r - a,) Qft? + ^AL, 



(33) 



dv ~ v r •**' *'■•.! '3 M 

aPr?j-i ) _ 3ffl«: 1 i ) i 

Consider first the case i^r. We may write 

where #$,/ is a homogeneous polynomial of degree y in w and t'. In fact, in 
order that this may be so, it is necessary and sufficient to satisfy the conditions 
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O. Tf(k) 



^±± +{(*-&) R?U = 0, 



3 -^r, i, 2 I /„ „ \ »») _ p(*-l) V R>r,i,1 _i_ / jO /0 \ D(t) _ /)(*-l) 

— g- — T V a r "-I) a r,i, 1 — r r,i,l) 5Tj r IPr Pi) "'r.i.X — Vr,«,l i 



o" T V« r a i/ a f, i/2k-Z — f f,t,!i-!l ^w ' ^ r Pi) IX r,i,ih~i — Vf,i,2*-2, 

(<* r a j)-^r, j, 2fc-l =: *>-, {, 2ft-l> (f'r Pi) -^r, i,2k—l == Qr,i,2k-1 > 

which are consistent on account of (33) and give 



T>(k) f f,i,it-l Vf. i, 2&-1 

3 p(fc-i) 



a.— a. 



*ft 



2fc-2 



" Vf,< ; 2t-1 

dv 
~d r -* (a r - ai ) 2 &.-& (&-&)» 



p«-l) 

. *■ r, i, 2fc — 



2J:-2 



9^ _gjfjI&-2 



p(fc) 



pOfc-1) 
■ *• r, i,2fc-X 



a, — a,. 



Vr, t",2fc-X 



3 p(*-i) 

(/ - r r,i,2t-Hl 

3w 



(a r — a,.) 5 

^ Vr,t,2ft-X+1 

3v 



+ .... + (-1) 



&-& ~ (&-&)» 



+ .... + (-1) 



\-i 



\-i 



(a r -a,.) x 
a\-i/i(fc-i) 

u V r, t, 2fc-l 



If we write 2Jc — X = j, we obtain the general formula 



»»> -"r 1 (-i) x a x nVAx _ " j- 1 (-i) x a x Q&? + x (34 x 

(j = 0,l,2, ....,2&-l), (» = 1,2,3,4; »=£r), 

which gives, upon substituting the values (32), the following two equivalent 
expressions : 

3X--1 /(«-l) 3X+1 A(k-1) pt\ A(.k-l)-\ 

\rr n> 3 w x-i du^dv du x ~ l dvy 



2&—j— 1 

*'**'' = x?o 4( { 



X^o 4(&-&)* +1 Kft-l) 

— * ( « r + a,) ° , ^, ,-+x-i - P a V' i +x - a, j ^ '+* , 
(; = 0,1, 2, ....,2fc-l), (* = 1,2,3,4; i^br). 



(35) 



31 



248 Wilczynski: On a Certain Completely Integrable 

For i = r, we find 

C = f ^ [SnVi } + • • • • +P£:&-i!dt*+ jg£r?+ • • • • H-G^-ii^L 

and the conditions (33) reduce to 
We may write 

'r, r — L-"r, r, + -Br, r, 1 + • • • • + -"r, r, 2fc] (a/*), 6 r W) j 

where i?<*>^ +1 is a homogeneous polynomial of order j -\- 1, subject to the 
compatible conditions 

dflffi.M-i _ p(fc-D _ i p 2r A (fc-D _ „ g^fr'H 

9« r ' r "' _ 4(a r r — 1)|_9** r, '~ 1 3^ J' 

U = l,2,3, ....,2*-l), 
where 

2r = -/3 f M 2 + 2«v-a f t; 2 . (36) 

If we multiply both members of the first equation by u, and those of the 
second by v, and add, making use of Euler's theorem on homogeneous functions, 
we find 

U + 1) BgU, = 4 iJ^ [2 fl,^ - (^^- + ^^)] > (37) 

(; = 1,2, ....,2*-l), 

to which we may add 

B r %i = 0. (37 a) 

The value of the arbitrary constant i?^ )0 is immaterial. The resulting ex- 
pression for y$ k) will be 

y m = e .,.+/»,. [B « s + . . . . + jjw u ] W)1iM 

+ e«'- + ""2' {B»U + • • • • + ^*U-i| 



4 



_ 2 / e *. + fl«. + c-,-..)..<» + tf,-fl.)»*« [ j R (*j j0+ .... +^ >2fc _ 1 ] 0lWj6iW , (38) 

and this will have the desired form, if and only if the limits are chosen in such 
a way that 

[■Rr.^O +••••+ Rr,i,2.k— l]o 4 <*>, M*> = ^> 



4 

[■"r, r, 2 + ••••.+ J*r, r, 2k 1 a r <*>, b r W = 2 H Tj {j . 

«=1 



(♦ = 1,2,3,4; izfzr), (39) 
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For every finite value of k this can certainly be done, and we shall then have 
the following formula for yf> : 

y m = ^.H-fl,. [R m 2fc + £ ( #«o. x + . . . . + #« 2fc _ l) j , 



4 

if we remember that, according to (37 a), 

K% i = o. 



Therefore, 

4#=:2fl r %, (i = l,2, ....,2*-l), A<*\ k = R<>l, th , 

i-l 

whence, making use of (35) and (37), 

1 * 2&-,-x / -i\x r a\ ^ (fc-x) 

+ j£ & (»,-a < )mft-l) [ i - (ft + ft) " + 2 "' g 4^= : 

ax-x -i(fc-x) ax+t ju-i) ax ^< fc - ] ) T 

0* = 1,2 2* — i). 

^ 2&_ 4fc(a r/ 3 r -l)^ 2fc - 2 - 

A second equivalent form for A%) may be found by using the second of the two 
equations (35). 

From the last formula we find at once 



^ = *»*l(«,ft-l)» ^ = 1 - 2 - 3 >■ («) 



Let us put 



^ x = £(a r - ai )* +1 K&-1) ' ^ x = £(a r - ai )^( a ;^-l) 



(41) 



Then we find 



r '' 4^(a r /? r — 1) L \ 3«* 3^ /J 

1 2*-,-x r ax j (&-x) ax-x j (fc-x> 

+ i & <- 1 ) > {<-'^ M + 2 ^^) g 4r = -'-"'-' at-'"" 

ax+x a (*-i) ax 4 (fc-i) T 

o r,, v thli ^+x o c r - i+h i f-icn 

6r>x ** du^dv ^^du^^dij ' (42) 

(; = 1,2,...., 2& — 1), 
and a second equivalent form corresponding to the second form of equation (35) . 
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Thus we may obtain a solution of (8) in the form 

y r = y°r + yPc + yf>c% +....+ yf^cl + m, 

where the further development of 31 would give powers of higher than the k-tli 
order in c , h being any finite integer, and each of the functions y® being of 
the form 

yf = e a '» +ft 'S A?). 

We may not, however, conclude in this fashion that a solution of (8) exists 
of the form 

2/r 0) + 2/r 1>c o + y? )( $> + .... ad infinitum, 

since we can not be certain that the limits of the integrals, as determined by 
equations (39), do not tend to become infinite with h. If this were the case, 
the convergence theorem which was proved in § 3 would not be applicable. 

In order to establish the existence of such solutions, we may proceed as 
follows: In equation (25) let us put a .= b = 0. We then find, making use 
of (30) and (30a), 

y m — e a r « + p r v P m ( UfV ) + ^ c a) e «. •+/».• f 

where 

1 i 

With this value of y^ we find further 

where 

<V> = -R<% , (s^r), 2 c"-" = 0, 
or 

y (B) _ e ^.+/» r . P ffl) (M)t; ) + j c (o 6 «.«+A»p(i) ( M>V ) _|_ | c« e a ' tt+f1 '", 

if we write 

4 

~(2) „(1, 1) I y r (l) r (.V> 

We may now prove by induction that 

y (») _ e a,.+Aipw ( W)V ) + 2 c o) e «.«+A»pM («, V ) + +2 c« e°' tt+ ^\ 

or 

y m — e a r u+e r v P m ( u ^ + £ | c co e «. «+/».. p<*H> ( M>v ), (36) 
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all of the quantities c®, being constants, and the functions P£ X) being poly- 
nomials of degree 2 X. 

We know that the series 

yr = e a r u + p TV+ | y (*, c J (37) 

h = l 

is absolutely and uniformly convergent for all finite values of u, v and c , and 
that for u = v = we shall have 

y r — *-> a — — a r> a — — "»•> 3—3- = a r p r . (do) 

From the general theory of linear partial differential equations, we know 
further that the uniquely determined solution of system (8), which satisfies 
the initial conditions (38), may be represented as a power-series $ (u, v,c ) 
in the three variables u, v, c , absolutely and uniformly convergent for all 
finite values of these variables. Moreover, for any system of finite values of 
u, v, c , there exists a convergent series D of positive constants which domi- 
nates the series of absolute values of the terms of $■ (u, v, c ). 

Since the series (37) and 5$ (u,v,c ) are both absolutely and uniformly 
convergent, and since both of these series represent, for all finite values of 
u, v, c , the same solution of (S), we shall obtain 5)5. (u, v, c ) from (37), if we 
substitute in (37) the expressions (36) for y ( r k \ expand the exponentials into 
power-series and perform the indicated multiplications. 

Now the quantities | e a ° u+ls ' v \ and | P£ X) (u, v) | are at most equal to the sum 
of the moduli of the terms of their respective expansions. Let us replace the 
several terms y ik) c%, of series (37), by the sums resulting from (36), and let 
us regard the individual terms of the sums thus obtained as terms of a new 
series 2. The series of the absolute values of the terms of 2 will be dominated 
by a series of positive constants, whose sum is at most equal to the sum of the 
series D of positive constants, which dominates the series of absolute values of 
the terms of ^ (u, v, c fl ) . Consequently, the series 2 itself and any series com- 
posed of terms selected from those of 2, provided that no term of 2 be taken 
more than once, will be absolutely and uniformly convergent. 

Therefore, the series 

CO 

y = e a r u + Pr v ^ p«) ( M)V ) c h f 
fc=0 

which is known to satisfy system (8) formally, is absolutely and uniformly 
convergent. If now we differentiate this series term by term, either once or 
twice, with respect to u or v, and compare the resulting series with the corre- 
sponding first- or second-order derivatives of ^ (u,v,c ), we recognize by a 
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repetition of the above argument that these series also are absolutely and 

uniformly convergent and therefore represent -^- , -?^- , etc. 

We finally obtain the following theorem : 

If u , v is a point for which the exponents a a , . . . . , a 4 are distinct, there 
exist four linearly independent solutions of system (S), of the form 

y r = e a,c«-«o)+ft-c»-»o) [l+| p(*) ( u — u , v—v ) c fc ] , (r = 1, 2, 3, 4), 

k=i 

where P*>, the coefficient of c$, is a polynomial of degree 2 k. 

In formulating this theorem we have made use of a terminology suggested 
by our discussion of the effect of the linear transformation (28a). In the 
vicinity of the point u Q , v , the coefficients of y in the two equations of (8) 
may be written in the form 

c (u — u ) + c Q u Q + c lt c (v — v ) + c v + c 2 , 

and we may speak of the four pairs of numbers a, /3 which satisfy the equations 

a 2 + 2/3 + c u + c, = 0, /3 2 + 2 a + c v + c 2 = 0, 

as the exponents for the point u , v . 

We see that, associated with a system of form (8), there is a set of poly- 
nomials. We have seen, moreover, how these polynomials may be calculated 
by means of certain recursion formulae. These formulae, however, are rather 
complicated, and we now proceed to study the polynomials P<*° (u, v) by a 
different method, which gives far simpler results. 

Let us substitute, in (8), 

yz=Yie au + fiv , 

where a and /? are any two numbers which satisfy the characteristic equations 

a» + 2 /3 + c x = 0, £» + 2 a + c 2 = 0. 
The system (8) will be transformed into the following system: 

y! nu + 2a.Y! u -\-2v! v + c q uyi = 0, y; vv + 2>? M + 2/3 >?„ + c v^ = 0. (40) 
Since (#) has a solution of the form 

fc=l 

(40) must have a solution of the form 

v! = 1 + 2 P (fc) (w, v) c fc , (41) 

where 

P<*> («*, v) = 2 ^ fc) («, v), (42) 



System of Linear Partial Differential Equations. 253 



A$ k) being a homogeneous polynomial of degree i. 

If we substitute the expression (41) for v\ into (40), we find 

Pit' + 2 aP<*> + 2 P<* > + ^P^" 1 ) = 0, 
P£> + 2P<*> + 2/3Pf> fuPC- 1 ) = 0, 

From (43) and (42) we find, further, 

2 ( ^Ft- + 2a^- + 2 ^- ) + M 2 ^f- 1 ) = 0, 
i=A or ou a v / ,-=i 



j (ft = 1,2, 3, ....)• (43) 



2ft— 2 



(44) 



2 f^|_ + 2^- + 2 0™±-) + v X A«-» = 0. 
If we equate to zero the terms of degree 2 k — 1, we obtain the equations 

9 « ^^ft _|_ O ^^-2ft _|_ .. J (ft-1) _ A n °^ 1 2fc _|_ O O t/^-2 fc I ., A (ft-1) _ A 

whence 

2(a/5-l)^ = ifl.^^, 2( a ^-l)^i > = |g )) J^-- 1) 1) , 
where 

and where a /3 — 1 is different from zero if, as we may assume, a is not a 
repeated root of the characteristic equations. These equations show that the 
Jacobian of A$ and q is equal to zero, so that A$ must be a function of q. 
We easily find 

^ = 4'(«A)'*» ' < t = 1 ' 2 ' 3 "-'-)' < 40 > 

a result which we had already obtained by our previous method. 

Since A^ k) = 0, the terms of degree zero and unity, in (44), require special 
attention. We find 

r) 2 A <*> 3/4 W 3 A W 3 2 /( {/c > 3 J(«) 3 /( (*) 

ow 2 9i* dv or om ov 

whence 

3 AW 32 JW 22 J (ft) 

2 (a/3-1) ^- = -|3 Vt- + -A- > 

1 9t« ' o% 2 ov 2 

a j (*) 32 4 (ft) 32 ^ (*) 

Ot) ou 2 ov 2 

Let us multiply the members of these equations by u and v respectively, 
and apply Euler's theorem on homogeneous functions. We find 

2 (a/* — 1) 
where fl denotes the operator 
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= { -p u+v) £?+ {u - av) &>=K q «£l? +q ^)> (47) 



which plays a fundamental role in this theory. In the same way we find, from 
(44), equating to zero the terms of the first degree, 

1 



AP = 



a (AP). 



(48) 



4(a/?-l) 

Let us now equate to zero the terms of degree 2k — i — Din (44), where 
1 < 2 k — i — 1 < 2 k — 1. A precisely similar treatment gives the following 
equation : 

2(2k-i) (*P-l)A2B- t = qAfell„ + Cl(A8L t+1 ). 
If we put 2 k — i = v, we find 

A ' =* 2v + (l(3-ir > (" = 1.2,3,. ...,2fc-l), (49) 

for this formula will also he valid for v = 1 or 2 if we equate to zero A!- k) 
whenever i becomes zero or negative. 

By means of formulae (40) and (49) all of the homogeneous polynomials 
Al k) may be calculated. We shall actually write down their values for 
k = 1 and k = 2, viz. : 






(50) 



and 



A?> = 

(-6a+ll/3 2 -9a 2 /3+5a* + 4a/3 3 -5^)u+(-6p+lla 2 -9ap 2 +&^ + 4:a»/3-5^)v 

32(a/J-l) 8 

A? = 

(-6 /3+5a 2 -±a/32 + 5l3*)u 2 + 2(2+ya/3-5a s -5p s +a 2 l3 2 )uv+(-6a+5p 2 -4 : a 2 p+5a*)vZ 
- 32(a/3-l)* 

A? = 

(2+3a/3-5/3*)u 3 +3(-ia+5P 2 -a 2 /3)u 2 v+3(-4:p+5a 2 -al3 2 )uv 2 +(2+3ap-5a 3 )v 3 

48(a / 8-l) 3 ' 

A?> = 

/3 2 u*-4:Pu 3 v + 2(a/3+2)u 2 v 2 -4:auv 3 +a 2 v 4 > _ q 2 

32(a/?-l) 2 32(a/J-l) 2 ' 

We observe in all of these cases that 

R<*) 

A (*) — D * 

' ~~ (a/3 — l) 3 *-"' 



(51) 



(52) 
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where B^ is a homogeneous polynomial of degree v in u and v, whose coefficients 
are integral rational functions of a, (J of degree 3 k — v. Making use of equation 
(49), it is easy to prove by induction that this law is true in general. 
If we make the substitution (52), our recursion formulae become 

B ®=Wk~V B ? ) = f,[ fl ( B i? 1 )+a^ 1) ], ( V = l,2, ....,2*-l), (53) 

where we must remember that B$> = B^X = 0. 

By repeated application of these formulas we find 



B (k) = 



2k— v 

2 ik ~"k\v(v + l) (2k— 1) + £ 2™-"~i +1 v{v + l) (2k— j) ' 



(v = l,2,....,2k-l), (54) 

where the symbol D. x indicates that the operation £1 is to be performed X times 
in succession. 

We may obtain another formula for 2?£ fc) from the equations 

B™ =l[q B<?-» + a (Bi%) ] , 

A V 



7?(*) — x 

r - 2 ~2(v—2) 



[qB£? + n(P£?>)], 



by eliminating the first terms occurring in the right members. In this process 
we must distinguish two cases, according as v is even or odd. We find 



X-l 



B& = 2 



g'ft^fcW 



B%U = 2 






X— 1 



iZo2*+ l (23,-1) (23,-3) (2A, — 2»— 1) 

or, to include both cases, 



(55) 



X-l 



2*11 (5£# +1 ) 



' 2 i+1 i/(y — 2) (v — 2i)' 

where 3, = ~ if r is even, and Jl = T" if r is odd. 

In the second equation of (55), let us put % = k. We find 

B m - v 1 ?& (g!l!) 

x>2 *- 1 - ^^-W (2A> — 1) (2fc — 3) .... (2& — 2% — 1) (* — t) ! 

Moreover, 

£1 (g*-«) = 2 (* — ») a* - * -1 *, + 4(k — i)(k — i — l)q> 
where 

32 



(56) 



fc — t— 2 



* 3J (57) 
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l 1= (p* — a )u+ (a 2 — P)v, 

l 3 = (! — £»)«» + 3 (£ 2 — a)u 2 v + 3 (a 2 — /?)mv 2 + (1 — a s ) v 3 



(58) 



Consequently, 



9 



fc — 2 



^ffi-^^-ra^iW+^M*)], (59) 

where 



fc— 1 Oi 

01 (fc) = ,? (2*— 1) (2ft— 3) .... (2ft— 2i— 1) (ft— *— 1) ! ' 

2 (ft) - 2 



(60) 



<=o (2ft— 1) (2ft— 3) (2ft— 2»— 1) (ft— t— 2) ! ' 

We may evaluate these sums and thereby very materially simplify equation 
(59) by the following method: If we equate to zero the terms of degree 2ft — 2 
in (44), a simple combination of the two resulting equations gives 

a--^-fl.-5?- 1 =(^- l )C ,, T#-- M -a#;- (61) 

In this equation let us substitute the values (53) and (59) for B$> and B^_ x . 
The left member of (61) becomes 

and the right member reduces to 

E = 4 "(ft — 1) ! ^""'^^-w^) + (ft — l)q h - 2 (vql — wg 2 )] . 
But 

( / 3 2 -a)^-(a 2 -/?)g„ = 2(a/?-l)( aM -,3<,), 

9 "l^ _& ll = 6(a/3_1)[_w3+(2a+ ^ 2)w2 ' ;_(2 ^ +a2)M ' ;2+ ' ;3] ' 

t;g|-tt^=4[-M 3 +(2a + ^)tt 2 v-(2/3 + a 2 )tti; 2 + i' 3 ], 
so that the equation L = R assumes the form 

q(au — ftv)[fli(ft) — (fe — l)l ] 

+ [_ tt 8 + (2a + /3«)««t;-(2/3 + a»)Mi;» + f;»][3 s (A;)— ** J = 0. 

The two cubic expressions which appear in the left member of this equation 
are linearly independent, since a @ — 1 is different from zero. Consequently 
we must have 
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*! (*) = 



so that 



(*-!)!■ 



M*) 



3 (A; — 2)!' 



B w _ g fc - 2 [3g^ + 2(fe-l)Z 3 ] 

^2*-l~ — 3 . 4 *( fc _ 1 )I 
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(62) 

(63) 



The explicit expression for P$L 2 , obtained by an extension of this method, is 
already rather complicated, viz. : 



B&- 2 = 



[18 S 2 (fc-2) -3 (a/3-1) \q s 



3 2 -4 k+1 (k — 2)\ 

+ 9\5ll+±(k—2)(aP—l)uv—8(aP—iyuv\q* 

+ 6\(Jc—2)(2l 1 l 8 —3l i )+3(k—2)(a{3—iyu z v 2 \q 

+ 4(k-2)(k-3)H], (64) 

where 

Z 4 = /?(— p u + vy+0L(u — av)\ (65) 

The transformation (28a) enables ns to simplify these formulae very con- 
siderably. Let us assume that such a transformation has been made so as to 
reduce c x and c 2 to zero. The characteristic equations become 

a 2 + 2/9 = /3 2 + 2a = 0, 

and one pair of solutions is a = (3 = 0. The operator fl reduces to 

3 2 3 2 

d« 2 or 

and it becomes an easy matter to calculate as many of the polynomials P (fc) (u, v) 
as may be desired, although, even with this simplification, the general expression 
seems to be rather complicated. We find 

P< fl > = 1, 
PW = — ~uv, 



p(») = -Atf» + ^(«' + « , )-|« , f' + i 8Wtp(tt , + » 8 ) 



P< 4 > = 



32 

7 



48 



« 3 f 3 , 



;MD 



2 6 -3 



(M s + v s )+ ■ u z v - 



2„,2 



^«i;(M 3 -)-t; 8 ) 



-J 1 QJ i 1)* 



(66) 



2 7 -3 
We shall close this discussion with one further remark. The expression 
(41) for r\ may be written 

(67) 



«o9 



^ = e *(afl-D + 2 Q<*> ( M , v) c\, 
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where the functions Q (fc) (u, v) are polynomials of degree 2 k — 1, differing 
from P (fc) (u, v) merely by the absence of the terms of degree 2 7c. Thus, by 
introduction of the first exponential term of (67), the polynomials P (&) have 
been deprived of their terms of highest degree. By means of equation (63) 
this process may be carried out one step further, depriving the polynomials 
<2 (fc) , in their turn, of their terms of highest degree. However, for the present, 
we shall refrain from any further developments in this direction. 

§5. Integral Equations Satisfied by the Solutions of (S). 

If Y is any solution of (S), it may also be regarded as a solution of the 
non-homogeneous system 

&«. + 2y„ + Cjy + c uF = 0, y vv + 2y u + c 2 y + c v Y = 0. 

The method of § 2 for integrating such a system is applicable, since the inte- 
grability conditions are all satisfied. Consequently we obtain the following 
relation : 



(J) Y--=c 2e*<«+e<*f T r_ 1 s— ivl \ ( fru + 2v)Y-u^\ du 



»£«-**:/ 4 ( ^,.-i) [{<-ft«+ 2 *i 

** a, 6 

+ [(2u-a i v)Y-v^}dv] 

4 

+ X a i e a * w+fiiV , 

i=l 

where a x , . . . . , a 4 are constants. The integrals on the right member of this 
equation are independent of the path of integration, as a consequence of the 
differential equations satisfied by Y, and have definite finite values since all of 
the solutions of (S) are integral transcendental functions of u and v. 

Moreover, if these integrals be denoted by l t , they have been determined, 
in accordance with the method of § 2, in such a way that, if we write 






we shall have 



just as though l lf . . . . , l i were constants. Consequently, the solution of the 
integral equation (J) will satisfy the same initial conditions for u = a, v = b 
as the function 

S a i e a * u+fi * v t 
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so that the equation of form (J) satisfied by y r will be obtained if we replace 
this sum by e artt+PrV . 

Let us change the notation for the variables of integration. Then, the 
integral equation satisfied by y r may be written as follows : 

(J r ) ¥ (u, v) = e a < u+(s * v 



+ C 



*f a, b 

Evidently the Liouville-Neumann method of successive substitution, as applied 
to this integral equation, will lead us back to the solution already obtained in § 3. 
It is easy to see, conversely, that if a t , /? 4 - are the four pairs of solutions 
(supposed distinct) of the equations 

a 2 + 2 § + e x — 0, /? 2 + 2 a + c 2 = 0, 
any continuous solution of ( J r ) , with continuous first and second derivatives, 
will be a solution of system (8), satisfying the given initial conditions for 
u = a, v = b, thus establishing, by a second method, the uniqueness of the 
corresponding solution of the integral equation (J r ). 
In fact, we find from (J r ) by direct differentiation 

Y uu + 2T V + Cl 7 = c _S r( —^—^ [(-& + 4 «- «,&!*) J + ft M |I 

— (l + 0|i»)|— — wJ^|-l. (68) 



But we have from (24) 



2 A, 



and $k t is the co-factor of a t @ t in the determinant A defined by (22). 
Consequently, 

X = 2 ai = 2 — = 2 a ' i ^ i — 4 

i=l Cl.iPi — 1 *=1 OLiPi — 1 <=i GtiPi — l " i=i <x< & — 1 " ? 

so that the right member of (68) reduces to — Co^Z. In the same way we 
may prove that I r must also satisfy the second equation of system (#). 

The integral equations (J) and (J r ) are simple examples of a type of 
integral equations which apparently has not yet been considered. The charac- 
teristic property of equations of this type is that, they contain open line inte- 
grals subject to the condition of being independent of the path. We shall leave 
the general discussion of integral equations of this kind for a future occasion. 
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"We shall close this section by pointing out a remarkable formula for the 
partial derivative, with respect to c , of any solution of system (8) . Let 

Co 

Then we see that yj is a solution of the non-homogeneous system 

/£,n >7w + 2i7, + ( c o« + c i) >? + uy = 0, 
rivv + 2Vu+ ( c o v + c i) >7 + vy = 0. 
The corresponding homogeneous system is (8). If we apply the method of § 2 

to the integration of system (S'), we obtain ~-^ expressed as a sum of integrals 

d c 

of exact differentials depending upon the solutions y lf ...,,y i of (8). Of 

course, ~-^ may also be regarded as a solution of an integral equation obtained 
o c 

from (J) by differentiation with respect to c . 
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